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Confocal microscopy of fluorescent labeled particles has been used to study the dynamical and 
structural properties of colloidal and granular matter in real space. Localization algorithms allow for 
a fully automatized determination of the three dimensional positions and translational motions of all 
constituent (spherical) particles in the observed volume. Though, particle deformations or rotational 
motions were hardly addressed. Here we present preparation and image processing techniques to also 
extract the deformation and rotational state of the particles. The deformation analysis is worked 
out for particles with a hollow sphere-like fluorescence intensity distribution that are more sensitive 
to small deformations. In the second case of rotations we utilize the angle dependence of the light 
absorption of the incorporated dye molecules to prepare optically anisotropic particles in analogy to 
polarized fluorescence after photobleaching. Rotations of these particles are expressed as intensity 
fluctuations in the confocal images. In contrast to existing methods our techniques do not reduce 
the quality of the actual particle localization. They can help understanding complex reorganization 
processes in arrested states of colloidal and granular materials during aging or under external stimuli 
such as shear or compression. 



I. INTRODUCTION 



The mechanical state of colloidal and granular sys- 
tems ranges from purely viscous suspensions to networks 
and concentrated systems exhibiting significant elasticity. 
Many parameter have an impact, including the choice of 
materials and sample compositions or external parame- 
ters like pressure and temperature. The dominant phe- 
nomena, such as yielding^ ^ or shear thinning^ are easily 
identified on a macroscopic level using diverse rheomet- 
ric method^^^. However, the microscopic understand- 
ing of these processes is often hindered due to a lack 
of appropriate experimental measurement techniques on 
a microscopic level. Diverse theoretical approaches and 
simulation method^^H^^ try to fill the gap but may show 
discrepancies due to the complexity of the system. 
Involved particle diameters range from 10 nm to 10 jiia 
for colloidal and from 1 /im to 1 mm for granular matter. 
These lengthscales favor optical investigation methods ei- 
ther in scattering geometry such as dynamic and static 
light scattering^^ or in direct microscopic imaging. Some 
of these imaging techniques are specialized to the extrac- 
tion of three-dimens iona l (3D) particle coordinates such 
as x-ray tomograph}^^^ which can be used for particle 
sizes down to few tens of micrometers. In this paper we 
deal with particle diameters down to 1 /im that can be 
investigated with confocal microscopy ^^^^ For loose gel 
structures as well as densely packed glassy systems 3D 
particle coordinates were determined with accuracies of 
about 0.05 /im^^. Phase transitions and coexistence re- 
gions were easily identified^^ Crystallization kinetics 
as well as its frustration due to irregularly shaped par- 
ticles or polydispersity were also investigated^^. Using 
confocal microscopes with high speed scanning^^ or mul- 
tiple beam illumination^ also fast dynamic proce sses, 
like e. g. colloidal flows under external sheap^^^, at- 



tracted attention. 

In contrast, colloidal gels features only slow particle re- 
arrangements during gravitational collapse^^ or under 
other external forces^^. In this and other kinetically 
arrested systems motions are suppressed because act- 
ing forces compensate each othei^^. If such systems are 
further compressed the loading force is heterogeneously 
transmitted through the particle contacts along so called 
force chains^^. The characteristics of this heterogeneous 
distribution of forces is currently under debate. It has 
been simulated^ and studied experimentally, e. g. via 
materials that show stress-induced birefringenc^^ in case 
of granular systems. However, smaller particles in the 
micrometer range and 3D samples are inaccessible with 
this method. Another approach to identify contact forces 
between touching particles uses soft spheres that deform 
under external force^^. In other words deformable par- 
ticles can be used to sense local forces. From a practical 
point of view, in particular hollow spheres are suitable, 
since their elastic constants can be tu ned by the shell 
thickness and choice of material d^^ l ^^ l When stressed, 
these hollow spheres deform to an ellipsoid-like object. 
The induced anisotropy and orientation of this ellipse 
can be used to quantify the strength and direction of the 
acting force. Here we present a method to automatically 
determine this deformation state of all hollow particles in 
the system together with their 3D coordinates from con- 
focal images. The analysis is based on an already exist- 
ing algorithm to locate fully labeled spheres by Crocker, 
Giiem and WeekPI. 

In case of non-balanced forces complex reorganization 
processes come into play that change the particle assem- 
bly. It is questionable whether these rearrangements in- 
volve rolling or sliding motions of the particles^^. Clearly 
frictional forces between particles in contact are essen- 
tial in this problem and already were investigated in 
AFM studies^^. Conventional confocal microscopy can- 
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not distinguish between rolling and sliding motions be- 
cause both processes manifest in simple displacements 
of the particles in the confocal images. Only additional 
information about their rotational motion will resolve 
this ambiguity and help to answer basic questions like: 
Can rotations be enforced for rough or irregularly shaped 
particles? How do particle-particle interactions influence 
rolling and sliding? 

Rotational motions have been studied experimentally via 
confocal microscopy for anisotropically shaped particles 
or groups of bound particles^^. However, such grouped 
particles might be hardly recognizable in dense suspen- 
sions. In addition, this approach immediately limits the 
choice of particle shapes. In particular spherical particles 
cannot be used. A possibility to resolve this shortcoming 
are particles with anisotropic optical properties. Spatial 
variations of fluorescent dye concentrations or incomplete 
metal coatings on homogeneously labeled particles^^ are 
directly reflected in the microscopic images and thus can 
be used to unambiguously determine the rotational state 
of the particle. However, at the same time such vari- 
ations lead to inevitable complications in the analysis 
of the particle coordinates and eventually an alteration 
of the interaction behavior. To circumvent this prob- 
lem, we induce an anisotropy by selectively bleaching 
dye molecules with specific orientations. This method 
was used in (polarized) fluorescence after photobleaching 
[(p)FRAPp for the investigation of the rotational dif- 
fusion of particles in solution^^ and the incorporation of 
dye molecules in membranes^^ It relies on the fact, 
that dye molecules preferentially absorb and emit light 
of specific polarization^. 

We introduce two new methods to get mechanical infor- 
mation of colloidal systems on a single particle level via 
confocal microscopy. Hollow spheres are used to sense 
local forces. A distinction between translational and ro- 
tational motions is possible with the help of optically 
anisotropic particles. The paper is organized as follows: 
After an introduction into the processing of confocal data 
and the localization of particles, the special case of hol- 
low spheres is discussed. The algorithm for extracting 
the state of deformation is presented and tested for a 
sample data set. In the last part we describe the method 
to prepare particles with an anisotropic distribution of 
dye molecule orientations. Rotations of these particles 
manifest in variations in the total fluorescence intensity 
if tested with linearly polarized laser light. 



II. MATERIALS AND EXPERIMENTAL SETUP 

For the deformation analysis non-labeled polystyrene 
(PS) template particles (1.97 fim in diameter and a poly- 
dispersity of 1.05) were coated with a 60 nm silica shelP^. 
Rhodami ne B w as covalently incorporated into the sil- 
ica matrid^^EH. The particles were dispersed without 
any further surface modification in a melt of modified 
PDMS (Laser liquid™with defined refractive index of 



n = 1.5780 ± 0.0002, Cargille Laboratories, USA) to 
match the refractive index of PS. 

In case of the rotational analysis silica particles with 
a diameter of 780 nm (polydispersity of 1.1) were syn- 
thesized with covalently bound Rhodamine E^^. Due to 
the stepwise growth of the particles the dye was incor- 
porated only in the core of the particles while the shell 
is unlabeled. The particles were dried in a observation 
cell made of glass and infiltrated with refractive index 
matching mixture of glycerol and water. The system was 
sealed with UV-curing glue to prevent evaporation. 

3D images were obtained with a home-made laser 
scanning confocal microscope (LSCM) system working 
in fluorescence mode [see Fig.flJ^^. The probing laser 
(Cobolt Samba™ 532 nm, 25 mW, Cobolt AB, Sweden) 
was focused with a 100 x oil immersion objective (Olym- 
pus, UPlanApo PH3, NA=1.35) into the sample and 
locally excited the contained dye molecules. A gal- 
vanometer based scanning unit (SCANcube® 7, SCAN- 
LAB AG, Germany) was used to scan the focus across 
the fixed sample in x,y-directions at a maximum speed 
of about two 2D-frames per second. The z-position of 
the focus was controlled by moving the objective with 
a piezo-positioning system (nanoX 200 S, ENV 40 SG 
nanoX, Piezosysteme Jena GmbH, Germany). The emit- 
ted fluorescent light [Rhodamin B (9-(2-Carboxyphenyl)- 
3,6-bis(diethylamino)xanthyliumchlorid) ] was collected 
by the same focusing objective and separated from 
the excitation laser light firstly by a dichroic mirror 
(Laser-beamsplitter z 532 RDC, 90% transmission above 
558 nm, AHF Analysetechnik AG, Germany) and directly 
before the detector by a laser clean-up filter (AFH Anal- 
ysetechnik AG, Germany). An avalanche photo diode 
(idlOO, idquantique, Switzerland) was used for detection 
at a maximum count rate of 17 MHz. Since single pho- 
tons were registered with high quantum efficiency only 
small input intensities of down to 50 nW were needed. 
For the polarization dependent bleaching of the sample 
a thin film polarizer (VIS 4 K, suppression ratio 1:4000, 
Linos Photonics GmbH & Co. KG, Germany) was placed 
between the dichroic mirror and the scanning unit. 



III. BASICS OF PARTICLE LOCALIZATION 

As the laser is scanned across the sample the de- 
tected fluorescence intensity maps the distribution of 
dye molecules inside the sample. If the dye molecules 
are incorporated into particles with a local distribution 
^dye(r) the detected fluorescence intensity Id at point 
r = (x, z) is proportional to 

/rf(r)~ I ffPSF(r-r') 5dye(r' - r„) j rfV (1) 

denotes the sum over all particles n at their respec- 
tive positions in the observation volume. The convo- 
lution of ^dye with the so called point spread functiorP^ 
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FIG. 1. Schematic setup of the home-made laser scanning 
confocal microscope. The polarizer between the dichroic mir- 
ror and the scanning unit is used only in the analysis of rota- 
tional motion. 



^fpsF accounts for the finite spatial resolution of the imag- 
ing process due to the diffraction limitation of light^^. 
The aim of any localization procedure is to extract the 
particle positions from the measured fluorescence in- 
tensity. Iterative deconvolution methods^^ can be used to 
solve equation (!]) for if the distributions ^dye and ^fpsF 
are known. While ^psF can be obtained theoretically^^ 
and in experiments^^, ^dye is accessible via reasonable 
assumptions. However, these algorithms suffer from long 
computation times. 

The problem is simplified for big and/or well separated 
particles for which ^^psf only marginally affects the ap- 
pearance of two neighboring particles compared to iso- 
lated ones. The convolution integral and the summation 
in equation ([T]) can be commuted with only small errors 
and the apparent distribution function 
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FIG. 2. Filtering procedure: b) In the initial confocal image 
the particles are represented by broad maxima at a low signal- 
to-noise ratio. After convolution of the intensity data with a 
background corrected cut Gaussian mask distribution, shown 
in a), the maxima are well separated and clearly identifiable 
in c). 



In our case the indexes refer to an orthogonal, equally 
spaced 3D grid. If the discrete mask distribution ^mask 
[Fig.[2]a)] is modeled to match the apparent distribution 
^app the 3D convolution intensity C is maximal in the 
center of the particle. Here the measured local fluores- 
cence intensity and the mask distribution overlap com- 
pletely. With increasing separation of mask and real par- 
ticle image the convolution intensity decreases. This pro- 
cedure not only sharpens the intensity maxima represent- 
ing the particles but also reduces data scatter by replac- 
ing Id with the locally averaged quantity C [Fig.[2]c)]. 
Now the particles can be easily localized by finding local 
maxima in C. An additional local background subtrac- 
tion further simplifies the localization process. This can 
be incorporated into ^mask by subtracting the reciprocal 
of its total sum: 



^mask — 5'mask 
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defines the appearance of a single particle in the confocal 
image. That is, the particles are represented by more 
or less pronounced and separated intensity maxima like 
in Fig.|2]b). In this case a time-expensive deconvolution 
procedure is unnecessary and generally is replaced by a 
more practical approacfP^. In an algorithm introduced 
by Crocker and Griei^ and extended by Weeks^^ a local 
filtering is used to improve the image quality and thus 
the actual particle localization. This filtering is done via 
an additional 3D-convolution of the experimental data 
set with a mask distribution ^mask- 
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-j,n-k)ld{ij,k) (3) 



In this step we switched from continuous coordinates r 
to discrete indexes z, j, k to account for the nominal spa- 
tial resolution of the experimental measurement. This 
set of indexes denote the position of the volume element 
at which the fluorescence intensity Id{hj^ k) is recorded. 



IV. RESULTS 
A. Localization of hollow particles 

In the following we will describe the peculiarities 
of localizing a particle with a hollow sphere intensity 
distribution^^. We found that by modeling ^mask to the 
apparent distribution ^app the above described filtering 
process can be used for particles with arbitrary spheri- 
cal symmetric profiles, i.e. also for hollow spheres. Also 
ellipsoid particles with moderate aspect ratios are acces- 
sible. The latter case will be important for the deforma- 
tion analysis in the following section. On the contrary 
anisotropically shaped particles without intrinsic mirror 
symmetries cannot be localized reliably. 

Although a fast implementation of the algorithm is 
already available in IDL^^ and Matlab computing lan- 
guages, the code for 3D-convolution was rewritten in IDL 
using 3D-fast Fourier transform calculation methodJ^. 



The original implementation can only cope with Gaus- 
sian intensity profiles since the 3D convolution is sep- 
arated into three one-dimensional calculations. So, in 
particular measured intensity profiles that have a maxi- 
mum at the rim of the particles cannot be handled. The 
full 3D calculation in our implementation of the convo- 
lution integral does not impose any restrictions to the 
mask profile ^mask but requires increased memory capac- 
ities and therefore is limited to smaller data sets. How- 
ever, the total computational afford does not change sig- 
nificantly if the whole 3D data set is split up into suffi- 
ciently small fractions in order to reduce memory occupa- 
tion. The coordinate data is recombined in the end. The 
part of the implementation for the final localization of 
the particles was not changed. The algorithm was tested 
for real data sets obtained with our confocal microscope 
described above. Hence, spatial resolution, background 
level and optical artifacts of the confocal images like e.g. 
in Fig.lslb) were representative. 



As an example of the general performance of the filter- 
ing procedure we present in Fig.plb) a 64x64 pixel^ sized 
subsection of a xy-slice of a complete 3D data set. The 
intensity mask ^mask in Fig.jsja) was modeled according 
to the appearance of the particles in the image. The con- 
volved image in Fig.[3]c) revealed well separated maxima 
at the hollow sphere positions that could be easily lo- 
calized. Two aspects are worth mentioning: Due to the 
mask normalization in equation Q and a lower cutoff at 
zero for the convolution intensity C the background was 
strongly reduced. Second, the maxima in the convolved 
image came out sharp with extensions comparable to the 
thickness of the particle shells in the original image. This 
becomes more obvious when looking at the intersections 
in Fig.[3]b) that do not cross all particles perfectly in the 
middle plane. Already small lateral (x,y) and/or verti- 
cal (z) deviations from the center lead to a significant 
reduction of convolution intensity. Moreover, the profiles 
in Fig.js]^) show less noise compared to the original data 
aiding the subsequent localization process. 

Any variation of the measured intensity distribution 
to the modeled mask profile ^mask gave rise to a broad- 
ened and distorted convolution signal. If the deviations, 
e. g. in particle size, were too large the maximal con- 
volution intensity was no longer located in the center of 
the particle but still had the shape of a hollow object. 
This is shown on the right side of Fig. [4]a) for computer 
simulated data. However, for relative variations in the 
extensions of test particle r and mask Vq below 25 % a 
localization was still possible. This was also true if the 
particles were no longer spherical but rather were de- 
formed to ellipsoids. Moreover, any broadening of the 
particle representations in the convolved image affected 
predominantly those directions where a deformation ac- 
tually was present. This directional broadening could 
be used to detect deformations of spherical particles and 
extract their orientations. 
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FIG. 3. xy-slice of the hollow sphere mask a), raw Id b) and 
convolved C c) 3D fluorescence intensity of hollow spheres. 
The color scale is linear and in arbitrary units. The line pro- 
files in the lower part of the images cross several particles. 
Even small deviations of the line from the center of the parti- 
cle lead to significant reduction of the convolution intensity for 
both, lateral (x,y) and vertical (z) directions. Hence, separa- 
tion distances are large in the filtered data set and positioning 
is easy and reliable. 

B. Deformation analysis of hollow particles 

After the localization of the hollow particles the de- 
formation of each particle was extracted individually in 
a two-step process. In the first step its orientation with 
respect to the laboratory coordinate frame was analyzed. 
Afterwards the actual extension of the particle was de- 
termined from the original confocal image. 

As already mentioned above deformed particles gave 
rise to an anisotropic convolution signal. This signal re- 
sembled the shape of a filled ellipsoid. The major axes 
of this ellipsoid corresponded to the main directions of 
the deformation and thus could be used to define the 
orientation of the deformed hollow particle. For an auto- 
mated determination of these major axes we introduced 
the intensity weighted inertia tensor 6^5. 

l^p^P p 

with a, b being one of the directions 7/, z and a, b as the 
intensity weighted mean positions. 

d = ^J-- -"^Ip-ap (6) 
l^p^p p 

The summations include all pixels p whose convolution 
intensity Ip exceed a pre-defined threshold value and lie 
in the region where the particle is localized. The defini- 
tion of Qab is analogous to the definition of the inertia 
tensor in classical mechanics with the mass replaced by 
the convolution intensity. By definition Qab is a symmet- 
ric matrix that can be transformed to its diagonal form 
O^^. The principal axes (m = 1,2,3) obtained from 
this mathematical procedure can be directly identified 
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FIG. 4. a) Deviations in the particle extension of the measured intensity r from the modeled mask To lead to a broadening of 
the convolved image expressed in the increase of Gnorm- Two xy-slices of the computer generated 3D particle images (top) and 
the convolution (bottom) are shown on the right. If the relative deviation are too large r/vo > 1.3 the convolution intensity is 
no longer maximal in the center, b) Deformed particles appear as filled ellipses in the convolved image. The eccentricity of these 
ellipses is quantified by Sxx/Oyy. This ratio is shown for varying deformations Ax/r in x-direction at a constant deformation 
in y-direction Ay/r. For matching mask and particle distributions r/ro = 1, Oxx/Oyy does not increase monotonously with 
increasing aspect ratio A = Ax/ Ay. This ambiguity can be resolved via a mismatch r/vo / 1. The arrows indicate parameter 
configuration for which real (top) and convolved (bottom) images are shown on the right. 



with the major axes of the ellipsoicF^ and thus defines 
the orientation of the deformed hollow particle. Yet, this 
procedure relies on the assumption that the particle de- 
formation is unambiguously transferred into the broad- 
ening of the convolution signal. As shown in Fig.[4]a) 
this was not always the case. 0^6 was calculated for the 
above discussed variation of the particle radii r at a fixed 
mask radius Tq. As the particles were still spherical in the 
convolved image [see rig ht in Fig.|4]a)] 6^6 was already 
diagonal with equal elements Qxx-, ^yy and Qzz- Nor- 
malized to the corresponding tensor element of the mask 
distribution ^mask, ©norm ^as minimal for r/ro = 1 and 
rose irrespective of whether the particle was smaller or 
larger than the mask distribution. Therefore, a particle 
that is compressed in one direction and stretched in the 
other can also lead to a spherical convolution signal for 
which the orientation cannot be extracted. 
This effect is demonstrated in Fig.[4]b). An initially 
spherical particle with radius r was varied regarding its 
x-deformation Ax/r at a fixed compression in y-direction 
Ay/r = 0.8. Here Ax and Ay denote the particle exten- 
sions in X and y direction, respectively. After the cal- 
culation of &ab the ratio of tensor elements Qxx/^yy 
gave a measure of the eccentricity of the convolution 
signal. If the eccentricity deviated from 1 the particle 
orientation could be determined. The mask distribu- 
tion was chosen to exactly fit the original particle size 
r/ro = 1. As expected Qxx/^yy equaled 1 for equal de- 
formation Ax/r = Ay/r = 0.8. But for a x-deformation 
of about Ax/r = 1.1 corresponding to an aspect ratio of 
A = Ax/ Ay = 1.5 the tensor elements evened out again. 
Hence no orientation could be extracted from the con- 
volved image although it was clearly visible in the initial 
image [see right in Fig.[4]b)]. 

In order to resolve this ambiguity the mask distribution 
properties were chosen to on purpose mismatch the par- 
ticle properties. This was tested in Fig.[4]b) for a size 
mismatch of r/ro = 1-14. The deformations Ay/r and 
variations in Ax/r were kept fixed for better compar- 



ison. For this particular set of parameters the eccen- 
tricity Qxx/^yy rose monotonically with increasing x- 
deformation allowing for an unambiguous determination 
of the particle orientation. As a rule of thumb the mask 
extension should be reduced to the smallest extension of 
the deformed particles. A variation of other mask dis- 
tribution properties like the shell thickness altered the 
slope of the curves in Fig.[4]b) and hence may be used to 
adjust the sensitivity in aspect ratios. Without further 
discussion we state that a orientation analysis is possible 
for moderate deformations. Although the actual num- 
bers depend on the confocal resolution, spatial sampling 
rate and particle properties, variations in the aspect ratio 
and deformation of up to 25 % can be handled without 
any restriction. 

Once the orientation of the deformed particle is known 
the actual deformation along the principal axes can be 
extracted. However, this cannot be done on basis of the 
diagonal tensor elements O^^ because these quantities 
are mathematically not fully decoupled. Instead the de- 
formation is measured from the original confocal images. 
The projection of the intensity distribution along each of 
the principal axes shows two distinct maxima resulting 
from the particle shell. The separation distance of these 
two maxima is a direct measure of the particle extension 
and also its deformation. 

The complete orientation and deformation analysis was 
tested for a real confocal data set of undeformed hollow 
spheres. More than 99 % of all particles, about 1000 in 
total, were located with a sub-pixel accuracy (< 5% of 
radius). After extraction of the principal axes of Qab the 
particles extensions were measured in the original un- 
filtered image giving the same mean value for all three 
directions [see Fig.jsja)] as it was expected for spherical 
particles. Beside the particle extensions the orientations 
of the principal axes was tested. Each of the three prin- 
cipal axes was assigned to the axes of the laboratory 
frame a with a maximal projection |u^ • a|. This pro- 
jection is a measure of the orient ational correlation of 
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FIG. 5. a) Normalized distribution of particle extensions for spherical particles. The extracted principal axes lead to identical 
mean particle extensions (dashed lines). As expected the orientational correlation with the outer coordinate system is flat and 
is centered around the theoretical value of 0.83 (dashed line in the inset), b) The analysis of the stretched and rotated data 
set reproduced the stretching factors of 1, 1.15 and 1.3. The main axes correlate with the rotated coordinate system axes (see 
inset). 



the principal axes and the laboratory frame. Due to the 
assignment of axes values below 1/a/3 ~ 0.58 do not oc- 
cur. Besides, the distribution of the projections is flat 
since there is no preferential direction for the principle 
axes. Indeed the experimental distribution in the inset in 
Fig.|5]a) was found to be rather flat and its average value 
agreed reasonably well the theoretical value of 0.83. 

In order to test the performance of the algorithm for a 
known deformation and rotation the same data set was 
stretched virtually along the x- and y-direction with ra- 
tios of Ax/r = 1.15 and A?//r = 1.3, respectively. Ad- 
ditionally, the data set was rotated to exclude a possi- 
ble influence of the choice of the laboratory system. In 
the statistical analysis the stretching ratios were repro- 
duced and the principal axes aligned with the axes of 
the rotated coordinates system giving a maximum in the 
distributions of the projections at 1 [see Fig.|5]b]. 

The only requirement for this analysis to work is a ho- 
mogeneous fluorescence intensity across the labeled parts 
of the particle since Qab is sensitive to intensity fluctua- 
tions. So, in principle also fully labeled particles can be 
used. However, here the variations in Qab for deformed 
particles are less pronounced handicapping the analysis. 

At this point we want to discuss some potential diffi- 
culties of the analysis that arise from the finite size of the 
PSF. If the extension of the laser focus largely exceeds 
the thickness of the labeled shell the measured intensity 
in the vicinity of two adjacent particles is enhanced. This 
directly affects the convolved image C as well as the cal- 
culation of Qab- An efficient solution to this problem is 
to clip the maximal measured fluorescence intensity by 
introducing an upper threshold. The situation is getting 
even worse if the shell thickness is comparable to the size 
of the PSF. In this case the shells of two adjacent parti- 
cles merge in the confocal images and the determination 
of the particle extension after the orientation analysis is 
hindered. This case should be avoided by an appropriate 
choice of particle properties and spatial resolutions. 
The confocal images are also distorted by the anisotropy 



of the PSF that is broader in z-direction compared to 
X and y-directions. As a consequence the particle ap- 
pears blurred in z-directions. We circumvent this prob- 
lem by blurring the whole 3D fluorescence intensity data 
in X- and y-directions via an additional convolution with 
a suitable Gaussian distribution. 



C. Preparation of orientation anisotropic particles 
via selective bleaching and theoretical polarization 
contrast 

In the following we discuss a method to quantify the 
rotational motion of a fluorescent particle. There are 
no demands on the shape of the particles but the dye 
molecules must be immobile, like e.g. rhodamine-B^ 
that is covalently bound inside the silica particles used in 
this study. We want to mention that dye molecules that 
are incorporated into polymeric particles by swelling with 
a CO- solvent cannot be regarded as immobile. Therefore 
such kind of particles are not suitable. The rotational 
analysis is based on the fact, that most dye molecules 
show an approxi mate d ipole characteristic during excita- 
tion and emissiorP^^. 

Assuming an isotropic orientation of chromophores 
before the bleaching process the absorption dipole 
moment^^ follows the normalized distribution 
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with /iA denoting the absolute absorption dipole moment 
and 1^ and (j) polar and azimuthal angles. In the presence 
of the electric field 



exp[z {ujt — k z)] 



(8) 



of a linearly polarized laser beam^the probability of one 
molecule to be optically excited is proportional to \ij,a ■ 
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FIG. 6. Bleaching geometry: The polarizations of bleach and 
probe beam enclose an angle a. The orientation of the ab- 
sorption dipole /J.A of the chromophore is homogeneously dis- 
tributed on the unit sphere parametrized by polar i9 and az- 
imuthal angles. The emission dipole /lie encloses a constant 
angle with /it a that depends on the chromophore. Therefore, 
all possible orientations of /itE lie on a cone around /k-a- 



EbI^'^. This relation is not only valid for the absorption 
of photons but also for the bleaching j)f chromophores 
if a single photon process is involvecP^. So, the laser 
beam bleaches predominantly those molecules with ori- 
entations parallel to the beam polarization. Moreover, if 
the bleaching is irreversible the number of chromophores 
decays exponentially resulting in an anisotropic distribu- 
tion 
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exp [-a(^B, Ma) cos(0)] (9) 



The proportionality constant a will be regarded as the 
bleaching strength and is a linear function in /i^, E'^ 
and the irradiation time t. 

During the subsequent confocal imaging with a probe 
laser beam 



Ep(a) = 



cos{a) 
sin(a) 




exi)[i{ujt-kz)] (10) 



enclosing an angle a to the polarization of the bleach 
beam Eb, the emitted fluorescence intensity is propor- 
tional to 

27r TT 



sm{^) d^d(t) (11) 

In the experimental realization the probe beam polariza- 
tion was adjusted by a polarizer that was placed between 
the dichroic mirror and the scanning unit in the confocal 
microscope ([T]). Since the excitation and collected flu- 
orescence light share the same optical path in this part 
of the setup the fluorescence intensity in (pTj ) is modi- 
fied by an additional factor |/L4e('^,^) • Ep| whereas /l^e 



denotes the direction of the emission dipole of the chro- 
mophore. /L^E and fiA are not necessarily parallel. In 
case of Rhodamine-B both textbftors enclose an angle of 
about 23 This restricts the number of possible orien- 
tations of /[Xe to a cone about /l^a (Fig.[6| which is taken 
into account by an additional integration: 



Z7T 7T 

Ii{a, a) ^ J J 



|Ax?(^,(/),a)-Ep(a) 







Z 7T 

J |[R(^,(/>)Axr*(23°,V^)]-Ep(a)|'# 



sin(^) di9d(l) 
with the combined rotation matrix 



(12) 



(R(^,(/>)),^ = 



cos(^) -sin(^) 

1 
sin(^) cos(i?) 

cos((/)) — sin((/)) 
sin(^) cos(0) 
1 



(13) 



The results of this calculation are illustrated in Fig.[7|(a) 
in dependency on the bleaching strength a and the rel- 
ative polarization angle a. The fluorescence intensity 
normalized to the initial state before bleaching could be 
fitted by 



rnorm 



(a, a) 



Ii{a, a) 



A/-abs[sin(a)]'' (14) 



m 



The variation of the fit parameters are depicted 
Fig.[7|(b). Whereas the exponent n varies moderately and 
is of minor interest, the fluorescence intensity contrast 
AI = /f^^™(a,0°) - /f^^™(a,90°) between bleached 
and non-bleached orientations needs to be optimized by 
choosing an appropriate bleaching strength a. If a is too 
small AI is small and the resolution in the particle ro- 
tation is bad. On the other side, if a is too large the 
fluorescence intensity in the bleached orientation is too 
small to localize the particle in the confocal image. 



D. Rotational analysis 

Instead of rotating particles we demonstrate the us- 
ability of this approach by changing the relative polariza- 
tion of bleaching and probe beam for immobilized par- 
ticles. Bleaching and confocal imaging was done with 
the same laser beam enabling for an easy control of the 
bleaching strength. After selecting a suitable sample 
area with about 500 particles a fraction of the particles 
were bleached to a level of about 30 % of the initial flu- 
orescence intensity equivalent to a total energy dose of 
0.16 /iJ/particle. After bleaching the laser intensity was 
reduced to prevent further significant bleaching during 
imaging of the sample. Fig. [S] compares confocal images 
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rotation angle a / (°) bleaching strength a 



FIG. 7. (a) Simulated ratio of the fluorescence intensity /b//nb of bleached and non-bleached particles as a function of the 
relative angle between the polarizations of bleach and probe beam for various belaching strengths a. The data was fitted by a 
biased abs[sm{a)]^ dependency, (b) With increasing bleaching strengths a the polarization contrast /b(90 °)//b(0 °) increased 
strongly (red shaded area) whereas the exponent n varied moderately. 



Unbleached a = deg a = 90 deg 




x-position / (|Lim) 

FIG. 8. Spatially filtered confocal images of a single layer of silica core-shell particles. Before bleaching the measured fluores- 
cence intensity per particle was constant. After bleaching the particles in the upper part of the image the intensity of these 
particles dropped down by a factor of about 5 compared to unbleached particles. After changing the relative angle between the 
polarizations of bleach and probe beam from a = 0° toa = 90° the detected intensity contrast to the unbleached particles 
decreased. 



of the unbleached particles alone and in relation to the 
bleached particles under relative polarization angles a of 
° and 90 ° . The differences in the fluorescence intensity 
are clearly visible. Beside to these two extreme values 
3D confocal images were taken for various relative polar- 
ization angles. From each image the position and total 
fluorescence intensity of each bleached particle was ex- 
tracted and averaged (Fig.[9|. The fluorescence contrast 
AI amounted to about 40 % and the functional depen- 
dency on a could be fitted to equation (14) when fixing 
the exponent to n = 2. 

Although the theoretical modeling disregards the strong 
focusing of bleach and probe laser beam by the objective 
the dependency on a is reasonably well described. The 
large statistical scatter in the relative intensities might 
be an indication of this discrepancy to the model as the 
divergence of the beam tends to blur the bleaching con- 
trast. 




T 

30 60 90 120 150 
relative angle / (°) 



FIG. 9. Rotation of the laser polarization leads to a variation 
of the fluorescence intensity that can be fitted by a sin^ be- 
havior. The intensity is normalized to the average intensity 
for all angles. 
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V. CONCLUSIONS 



We showed that fluorescent confocal microscopy can- 
not only be used to obtain 3D coordinates of a particle 
but at the same time can deliver information about its 
deformation and rotation. This can be done on a single 
particle level in a fully automated procedure. Although 
there are no specific demands on the particles in partic- 
ular hollow spheres are a good choice for the investiga- 
tion of deformations in colloidal and granular systems. If 
the shell is thin compared to the used excitation wave- 
length {d < A/10) there is no need for exact matching of 
the refractive indexes of particle and solvent. Hence, the 
choice of compatible materials is considerably increased. 



On the other side rotational and frictional sliding mo- 
tions of particles have a strong impact on reorganization 
effects in arrested colloids. The described method can be 
used to extract this information, even for dense systems 
of spherical particles. 
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